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Abstract 



Using data collected with the CLEO II detector, we have observed the isospin- 
violating decay D* + Dfir . The decay rate for this mode, relative to the 
dominant radiative decay, is found to be T(D* + — ► Dfir°)/T(D* + — ► Dfj) = 
0.062t°;° 2 °± 0.022. 

PACS numbers: 13.25.Fc, 14.40. Lb 
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Since the discovery of the charmed strange vector meson D* + , only the radiative decay 
mode D* + — > Df~{ has been observed [1,2]. Due to the small mass difference between the 
D* + and the Df , the only kinematically allowed strong decay is D* + — > Dfir . This decay 
violates strong isospin conservation, since the initial state has 1 = and the final state has 
1=1. However, isospin is not an exact symmetry, so this decay is not completely forbidden. 
For example, the decay ip' — > J /if) 7r°, which has been observed by several experiments [1], 
violates isospin conservation. 

Cho and Wise [3] have made a prediction of the decay rate for D* + — > Dfir . Using Chiral 
Perturbation Theory, they describe the decay as an isospin-conserving decay involving a vir- 
tual ?7, D* + — > Dfrj; the r\ couples through its ss component, so this decay is not suppressed 
by the Okubo-Zweig-Iizuka rule. This is followed by the r\ mixing into a tt°. This isospin- 
violating mixing vanishes in the limit of equal u and d quark masses. The decay amplitude 
is proportional to the light quark masses in the combination (m^ — m u )/[m s — (m^ + m„)/2], 
which is ~ 0.02-0.03 [1]. They also conclude that the radiative decay rate is suppressed (rela- 
tive to D*° — > -D°7) because of the partial cancellation of the magnetic moments of the charm 
and strange quarks (this type of cancellation also accounts for the small radiative decay rate 
of the D* + [4]). Thus Cho and Wise relate the rate for D* + -> D+7 to that for D* + -> £>+7 
in order to estimate the ratio of partial widths: R = T(D* S + — > Df tt°) /T(D* s + — > Df^f) = 
~ 0.01—0.10. Unfortunately, there are corrections to the prediction which might be negligi- 
ble, but are presently uncalculable. There is also an electromagnetic amplitude for the decay 
to D+7T , but it is expected [3] to be smaller than the strong amplitude by a factor of order 

We have searched for D* + — > Dfir using the decay chain Df — > </>7r + , <f> — > K + K~. At 
the same time, we observe the radiative decay to normalize the hadronic decay rate. Since 
the same charged-track selection criteria are used for both decay modes, only the relative 
efficiencies for finding a single photon or reconstructing a tt° are needed. 

The data used in this analysis were collected with the CLEO II detector at CESR. The 
detector consists of a charged particle tracking system surrounded by time-of-flight (TOF) 
scintillation counters. These are in turn surrounded by an electromagnetic calorimeter which 
consists of 7800 thallium-doped Csl crystals. The inner detector is immersed in a 1.5 T 
solenoidal magnetic field generated by a superconducting coil. Finally, the magnet coil is 
surrounded by iron flux return and muon counters. Charged particle identification is provided 
by specific ionization measurements in the main drift chamber and by TOF measurements. 
A detailed description of the detector can be found elsewhere [5]. 

The data were taken at center-of-mass energies equal to the masses of the T(35') and 
T(4S'), and in the continuum above and below the T(4S'). The total integrated luminosity 
is 3.75 fb _1 . Events were required to have a minimum of three charged tracks, and energy 
in the calorimeter greater than 15% of the center-of-mass energy. Charged tracks were 
initially required to pass a loose particle identification consistency. We required that the 
specific ionization measurement be within three standard deviations of that expected for the 
hypothesis in question, either kaon or pion. 

Only energy clusters in the barrel calorimeter with | cos#| < 0.71 (where 6 is the polar 
angle with respect to the beamline) which were not matched to tracks were used as photons. 
They were required to have a minimum of energy of 30 MeV, and to pass a lateral shape 
cut to help eliminate energy from hadronic interactions. Single photons used to reconstruct 
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the radiative decay were required to have energy greater than 50 MeV. Pairs of photons 
were used to reconstruct 7r°'s. The invariant mass of the two photons was required to be 
within 2.5 standard deviations of the 7r° mass; this cut takes into account the asymmetric 7r° 
lineshape and the small momentum dependence of the mass resolution. The 7r° candidates 
were kinematically fit to the 7r° mass to improve momentum resolution. The decay angle, 
# 7 , is defined as the angle between the direction of one of the photons in the 7r° rest-frame 
and the 7r° direction in the lab-frame. We required |cos# 7 | < 0.75, since the background 
peaks near | cos # 7 | = 1 while the signal is flat. 

We began the reconstruction by taking pairs of oppositely charged tracks, consistent 
with being kaons, and calculating the invariant mass. Those pairs whose invariant mass was 
within ±9 MeV/c 2 of the <f> mass [1] were accepted as <f> candidates. Each remaining charged 
track, consistent with being a pion, was combined with the <f> to make a Df candidate. 
The Df candidates were required to pass two angle cuts. To reduce background from slow 
pions, we required cos 6 V > —0.9, where 6 V is the decay angle of the tt + (the angle between 
the pion's direction in the Df rest frame, and the Df } s direction in the lab frame). The 
signal distribution is flat, while the background peaks near cos 6 V = —1. Second, because 
the <f> is polarized in the helicity-zero state, the kaons must have a helicity angle distribution 
proportional to cos 2 6k, where Ok is the angle between the kaon and the Df , both measured 
in the <f> rest-frame. We required | cos 6k\ > 0.35. 

Next, we imposed a more restrictive particle identification cut on the three-track com- 
bination. A particle ID % 2 was calculated using the specific ionization measurements for 
each of the three tracks and the TOF measurements for each track which had good TOF 
information. We required that the % 2 probability be at least 0.1. Finally, Df candidates 
had to have a mass within two standard deviations (±16 MeV/c 2 ) of the Df mass [1]. 

To reconstruct the radiative decay mode, the Df candidates were combined with each 
photon in the event which had an energy of at least 50 MeV. The scaled momentum, x, 
of each D* + candidate was calculated as x = pjpu, where p 2 M = Eq — M^ t+} and E is 
the beam energy; we required x > 0.6. The mass difference, AAf 7 = M(Df^f) — M(Df) 
was calculated and histogrammed [6]. The resulting distribution was fit using a Gaussian 
modified with an enhanced low-energy tail for the signal and a third-order polynomial for 
the background. We find 944 ± 57 signal events (statistical error only). 

Next, each Df candidate was combined with each 7r° with momentum of at least 
250 MeV/c. An x cut of x > 0.6 was again applied to each D* + candidate. In Fig. 1, 
we show the mass difference, AM V = M(D+tt°) — M(Df) for the remaining D* + candidates. 
The data were fit using a Gaussian for the signal and a square-root function that goes to 
zero at threshold to represent the background. The r.m.s. width of the Gaussian was fixed at 
a = 1.2 MeV/c 2 , as determined by the Monte Carlo. The mean was fixed at 144.22 MeV/c 2 , 
as previously measured by CLEO [7]. We find 14. 7^0 signal events (statistical errors only). 

We also studied the Df and 7r° sidebands. We used </>7r + combinations with masses 
between 1904 and 1936 MeV/c 2 , and between 2004 and 2036 MeV/c 2 , and combined them 
with 7r°'s. We also selected 77 combinations that were between 2.75 and 7.75 standard 
deviations away from the 7r° mass; this corresponds approximately to an invariant mass 
between 88 and 118 MeV/c 2 , or between 145 and 165 MeV/c 2 ; these were combined with 
the Df candidates. These two sets of sideband combinations produce the AM V distribution 
shown as a dashed histogram in Fig. 1. The entries in this histogram have been scaled by a 
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factor of 0.5, to account for the fact that the sidebands are twice as wide as the signal band. 
If we fit this sideband histogram the same way as the signal band events, the area of the 
Gaussian is — 1.0l 2 ; 4 events, consistent with zero. 
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FIG. 1. The distribution of the mass difference, AM,, for the isospin- violating decay 
D* + — >■ _D+7T°. The points are the data after all cuts, the solid line is the fit to the data, and 
the dashed line is an estimate of the backgrounds, as described in the text. 

As a check, we have also analyzed events in the Cabibbo-suppressed decay chain D* + — > 
D + tt° } D + — > <f>TT + . We used the same cuts as before, but required the </>7r + invariant mass 
to lie in the D + signal region from 1853.3 to 1885.3 MeV/c 2 . Using a Gaussian with mean 
fixed at 140.64 MeV/c 2 , as measured by CLEO [8], we find 28. 6l 5 ; 4 signal events (statistical 
errors only). Using our measured luminosity, and published cross sections [9] and branching 
ratios [1], we expect 23.7 ± 4.2 events. 

In order to confirm that this signal is from D* + — > Dfir rather than D* + — ¥ -Dj"77 we 
have relaxed the cut on the two photon invariant mass used in the 7r° selection. We then 
selected events with 141.22 MeV/c 2 < AM, < 147.22 MeV/c 2 , and studied the 77 invariant 
mass distribution. Fitting this distribution yields a tt° signal of 16.lt^o events, consistent 
with our previous result. Similarly, when we cut on AM T and fit the <f>Tr + invariant mass 
distribution we find 13.9^4^ signal events; this is shown in Fig. 2. 



5 



The background appears to be dominated by random combinations, rather than feed- 
through from some other physics channel. For example, two other conceivable sources of 
background, real D* + — ¥ Df"f events with an extra soft photon faking a 7r°, and misidentihed 
D* + — > D + tt° events are both negligible. We have generated D* + — > Dfj Monte Carlo events 
and analyzed them with the reconstruction program. The Monte Carlo sample is 50% larger 
than our actual data sample. Only 3 events with 135 MeV/c 2 < AM V < 160 MeV/c 2 are 
found; they are all outside the signal region 141 MeV/c 2 < AM T < 147 MeV/c 2 . 



8 



6 i 



CM 



> 



V) 

c 

o 
> 

LU 



4 -II 



2 J 



0040795-002 



T 




1900 



1940 



1980 



2020 



M (ip7r + ) ( MeV / c ) 



FIG. 2. Fhe M(^>7r + ) distribution for events in the AM V signal region. Fhe points are the data 
after all cuts; the solid line is the fit to the data, using a Gaussian of fixed mean and width for the 
signal, and a second-order polynomial for the background. 



Similarly, background from misidentifed D* + decays is not a problem. The particu- 
lar decay chain considered is: D* + —* D + tt° } D + — > K~tt + tt + } where one of the pions is 
misidentihed as a kaon, so that the three charged tracks reconstruct near the Df mass. 
These events are not a problem because it is almost impossible for such a fake kaon to make 
a </>, and then for the fake <f> plus the real pion to make a Df. This was tested by taking 
each Df candidate and changing the particle identification of the kaon with the same sign 
as the pion, from kaon to pion. We recalculated the invariant mass of the the three tracks 
and found that none of the calculated masses were greater than 1845 MeV/c 2 . Studies of 
the data and Monte Carlo events confirmed that such D* + events could not produce this 
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peak. Other similar decays of charmed mesons, which are partially reconstructed, and in- 
clude misidentified particles, can contribute to the signal histogram. Monte Carlo studies 
showed, however, that such events do not form a peak. 

We have evaluated the statistical significance of D* + — > Dfir events in two ways. First, 
we refit the signal histogram (Fig. 1), constraining the area of the Gaussian to zero. The 
change in likelihood from the original fit is equivalent to 5.0 standard deviations. 

Second, we counted the number of events with 142 MeV/c 2 < AM V < 146 MeV/c 2 in 
both the signal histogram and the sideband histogram. This yields 16 signal events and 
5 sideband events. Taking into account that the sideband width is twice the signal band 
width, the binomial probability of getting 16 (or more) signal events out of 21 total events is 
7.3 X 10~ 5 , equivalent to 3.9 standard deviations. Thus we consider the observation to have 
at least 3.9 standard deviation significance. 

Using the CLEO Monte Carlo program, we found that the ratio of efficiencies for recon- 
structing the hadronic and radiative events is 0.25 ± 0.03. Using this ratio and the numbers 
of events reconstructed, we find that the ratio of partial widths is R = 0.062+Q;Q 2 g, where 
these statistical errors are dominated by the error on the number of Dfir events. 

We estimate that the systematic error on R is 35%. This is dominated by variations in 
i?o when we vary our cuts on the 7r° and/or photons. Thus our measurement of the ratio of 
partial widths is [10]: 



t(d:+ -> d+ 



If we assume that these two branching fractions sum to one, the individual branching frac- 
tions are B(D*+ -> Dfir ) = 0.058tffi ± 0.020 and B(D*+ -> D+ 7 ) = 0.942t°;°i8 ± 0-020. 

The observation of this decay mode implies that the D* + must have natural spin-parity 
(0 + , 1~, 2 + . . .), since conservation of parity and angular momentum forbid the decay of a 
particle with unnatural spin-parity to two pseudoscalars. The radiative decay rules out ± ; 
the most likely spin-parity is J p = 1~, the same as the D*° and D* + [1]. 

Using the D* + — > Dfir events, we can also make a measurement of the mass difference 
M D *+ —M D + , and set an upper limit on the width of the D* + . We refit the AM V distribution, 
allowing the mean of the signal Gaussian to float. With the width of the Gaussian fixed, the 
mean is fit to AM, = 143.76±0.39 MeV/c 2 (statistical error only). Fitting the D*+ -> D+vr , 
D + — > <f>TT + events mentioned above, we find Mz>*+ — Mz>+ = 140.31 ± 0.26 MeV/c 2 , in good 
agreement with the previous CLEO measurement of 140.64±0.10 MeV/c 2 [8]. We include the 
0.33 MeV/c 2 difference in the systematic error of the AM, measurement. It has previously 
been estimated that the uncertainty in the crystal energy calibration introduces a systematic 
error of 0.04 MeV/c 2 in this type of measurement [8]. Changing the cuts used to select the 
events introduces variations of 0.22 MeV/c 2 ; other effects, such as varying the background 
function used for the fit, produce much smaller variations. Thus we estimate the systematic 
error to be 0.40 MeV/c 2 . Therefore we measure M D *+-M D + = 143.76 ± 0.39 ± 0.40 MeV/c 2 , 
in excellent agreement with the previous CLEO value 144.22 ± 0.47 ± 0.37 MeV/c 2 from the 
radiative mode [7]. The two measurements are statistically independent, and have almost 
completely independent systematic errors. Averaging the two, we find M D *+ — M D + = 
143.97 ± 0.41 MeV/c 2 . This may be compared with the Particle Data Group's fit value of 
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141.6 ± 1.8 MeV/c 2 , or their average value of 142.4 ± 1.7 MeV/c 2 [1]; the PDG's values do 
not include the earlier CLEO measurement. 

The width of the signal is consistent with being entirely due to detector resolution. 
The Monte Carlo calculation predicts a value of a = 1.19 ± 0.07 MeV/c 2 for the signal 
Gaussian (statistical error only). The measured r.m.s. width of the signal Gaussian is a = 
1. 061^28 MeV/c 2 . Assuming a 10% systematic error on the Monte Carlo prediction for a, 
and fitting the signal with a p-w&ve Breit-Wigner convoluted with a Gaussian, we can set 
a 90% confidence level upper limit T(D* S + ) < 1.9 MeV/c 2 . This can be compared with the 
best existing limit, from ARGUS [11], T(D*+) < 4.5 MeV/c 2 . 

In conclusion, we have detected the isospin- violating decay D* + — > D+7T , and find the 
ratio of partial widths, R = 0.062+Q; 2 g ± 0.022, which confirms a recent prediction [3,12]. 
We determine the branching fractions for the Df^f and Dfir modes assuming that any 
other decay modes of the D* + are negligible. In addition, we present a new measurement 
of M D *+ — M D + } and we obtain an improved upper limit on the width of the D* + . The 
observation of this decay mode implies that the D* + has natural spin-parity, most likely 1~. 
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Sloan Foundation. 
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